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Abstract 

We have considered the propagation of gravitational waves (GW) in de Sitter space- 
time and how a non-zero value of the cosmological constant might affect their detection 
in pulsar timing arrays (PTA). If A ^ waves are anharmonic in Friedmann- Robertson- 
Walker coordinates and although the amount of this effect is very small it gives noticeable 
q \ effects for GW originating in cxtragalactic sources such as spiraling black hole binaries. 

O h The results indicate that the timing residuals induced by gravitational waves from such 

| \ sources in PTA would show a peculiar angular dependence with a marked enhancement 

around a particular value of the angle subtended by the source and the pulsars, depending 
mainly on the actual value of the cosmological constant and the distance to the source. 
The position of the peak could represent a gauge of the value of A. The enhancement 
that the new effect brings about could facilitate the first direct detection of gravitational 
waves while representing a local measurement of A. 
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1 Introduction 



Pulsar timing arrays (PTA) are one of the most promising candidates to offer the first direct 
detection of gravitational waves. They have been collecting data already for almost a decade 
and they are expected to obtain signals in the next years. The idea behind PTA is to detect 
the correlated disruption of the periods measured for a significant number of pulsars due to 
the passing of a gravitational wave through the system [IJ El [3j |4j. The frequency range 
sensitive to this method is 10 _9 s~ 1 < w < 10 _7 s _1 [I], and the timing residual is expected to 
follow a power law 0[5]. A key problem in making predictions for these signals is modeling 
in a realistic way the wave functions produced in the different sources, in particular the value 
of the amplitude of the metric perturbation h is a free parameter in principle. Some bounds 
in the range of 10~ 17 < h < 10~ 15 have been set already [5]. 

If A / gravitational waves (GW) propagate in a de Sitter space-time not in flat 
Minkowskian space-time. The general practice is simply to account for the expansion of 
the universe by using a redshifted frequency according to the distance of the source [3J. In 
this work we go beyond this exceedingly simple approximation and use an approximate so- 
lution of the GW equation in de Sitter previously derived [6] and see that the conclusions 
change. 

We assume that A is somehow an intrinsic property of space-time rather than an effective 
description valid at extremely large scales. If so, it is expected to be present at virtually 
all scales, with the possible exception of gravitationally bound objects such as galaxies or 
local groups of galaxies. If A is a fundamental constant of nature surely there should be a 
way of determining locally its value. By 'locally' here we mean at redshifts z « 1. This 
question has been addressed in [8] with varying conclusions. We will see that GW may 
open a nice window to realize this program. In fact, our results suggests that the currently 
observed non-zero value of A may actually facilitate the first direct detection of GW under 
some circumstances. 

This paper is organized as follows. In Section 2 the wave functions used are presented, 
the way in which the timing residuals are calculated is defined and a brief explanation of 
the coordinate systems involved is included. Section 3 is devoted to present our numerical 
analysis. In Section 4 we discuss the possibility of using this method to get some results on 
the value of the cosmological constant. Finally we sum up our conclusions in Section 5. 



2 Gravitational waves and timing residuals with A ^ 

In Minkowski space-time, gravitational waves obey the simple wave equation D/i = 0. It 
is possible to show [6] that in de Sitter space-time with A ^ and within the linearized 
approximation one can find solutions of the linearized Einstein equations in the traceless 
Lorenz gauge (TT gauge [9]) which obey the same equation of motion 

n^f = o. (i) 

Spherical massless waves are solution of this equation away from the source 

hH s = - (E^ cos[w(t - r)] + Dp, sm[w(t - r)\) . (2) 

However, as shown in [6], this simple linearized solution only holds in a specific set of co- 
ordinates, the Schwarzschild-de Sitter (SdS) coordinates. This is easily seen by considering 
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a linearized background solution (rather than wave-like solutions) and realizing that their 
unique static solution is the (linearized) Schwarzschild-de Sitter metric [B]. 

Although providing a perfectly valid solution for gravitational waves, SdS coordinates 
are not adequate to make observational predictions. The proper isotropic and homogeneous 
coordinates are the Friedmann-Robertson- Walker (FRW) onea3 and the solution ([2]) in such 
coordinates, neglecting 0(A) and higher, reads 



h FRW _ E /J.f 

" R 



cos 




w(T -R) + - TR 



sin 



w (T-R) + w x lU^-TR 



(3) 



where R is the usual radial FRW comoving coordinate and T is cosmological time. Note 
that the linearization process that has been used makes sense as long as AT 2 , AR 2 « 1 
and also that in the TT gauge the only spatial components of the metric that are different 
from zero are the X, Y entries of the polarization tensors E^, D^ u . Although some temporal 
components of E^v and are also non-zero in these coordinates, they are several orders 
of magnitude smaller than the spatial ones and therefore will not be relevant for the present 
study. 

We note that the phase velocity of propagation of the GW in such coordinates is not 
v p = 1 but v p ~ 1 — J±T + 0(A) [6j. On the other hand, with respect to the ruler distance 
travelled (computed with gij) the velocity is still 1 (up to terms in A of higher order to those 
considered) . 

Consider the set up depicted in Figure CD describing the relative situation of a GW source 
(possibly a very massive black hole binary), the Earth and a nearby pulsar. The timing 

Source 




Pulsar 



Earth 



Figure 1: Relative coordinates of the GW source (R=0), the Earth (located at Z = Ze) with 
respect to the GW source and the pulsar located at a coordinates P = (Px,Py,Pz) referred 
to the source. The Z direction is chosen to be defined by the source-Earth axis. Angles a 
and (3 are the polar and azimuthal angles of the pulsar with respect to this axis. 



1 Note that the FRW metric cannot be approximated to obey any linearized Einstein equation, see [6] for 
a detailed discussion. 
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residual [10] induced by (J3]) will be given by 



T r L 

H(T E ,L,a,f3,Z E ,w,e,A) = -— h % n> \ dx hf™ (T E + -x, P + L(l + x)n) (4) 

along the null geodesic from the pulsar to the earth, where we assum^l e ~ \E-ij\ ~ \Dij\, 
i,j = X, Y and the unit vector h is given by (— sin a cos /3, — sin a sin /3, cos a), Z# is the 
distance from the Earth to the source, L the distance to the pulsar and T E the time of arrival 
of the wave to the local system. In deriving the previous timing residual we have neglected 
the (non-zero) time components of E^^D^y that, as previously indicated, are several orders 
of magnitude smaller. The speed of light has been restored. We have assumed that from the 
pulsar to the Earth the electromagnetic signal follows the trajectory given by the line of sight 
R{x) = P + L(l + x)h. Since we assume that within the Galaxy A = 0, L is also the ruler 
distance. Explicitly 

R(x) = P + L(l + x)h = (—xL sin a cos /3, — xL sin a sin /?, Z E + xL cos a) (5) 
or in modulus 

R(x) = J Z E + 2xLZ E cos a + x 2 L 2 ~ Z E + xL cos a, (6) 

since we are considering L « Z E . This approximation does not affect in any significant way 
the results below. We do not consider here the known contribution to the timing residual H 
from the Earth's peculiar motion either. The integral is of course independent of the angle 
/3 for any single pulsar but it will depend on the relative angles when several pulsars are 
averaged. 

Let us consider the arguments of the trigonometric functions in (3) and define 
Q(x, T E , L, a, f3, Z E ,w, A) 



Then 

H(T E , L, a, /3, Z E ,w, e, A) = (sin 2 a cos 2 (3 + 2 sin a sin /3 cos 2 /3 — sin 2 a sin 2 0) 

2 c 

f dx <7 I X j T ( 1 + \^ Te + ~ x) J (cos + sin 0) • 

J_i (Z E + xL cos a) V V o c / 

(8) 

At this point one should ask whether the observationally preferred exceedingly small value 
of the cosmological constant [TT] affects the timing residuals from a pulsar at all. To answer 
this question we take reasonable values of the parameters both for the GW and one pulsar 
location and plot a snapshot of the resulting timing residuals as a function of the angle a for 
the time of arrival of the wave to the local system, T E . The comparison is shown in Figure 
[2j The figure speaks by itself and it strongly suggests that the angular dependency of the 
timing residual is somehow influenced by the value of the cosmological constant, in spite of 
its small value. Another feature that catches the eye immediately is an enhancement of the 
signal for a specific small angle a (corresponding generally to a source of low galactic latitude, 

2 This approximation is unessential and can be easily removed. 
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Figure 2: On the left the raw timing residual for A 
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-35, 



as a function of the angle 



a subtended by the source and the measured pulsar as seen from the observer. The figures 
are symmetrical for ir < a < 2ir. On the right the same timing residual for A = 0. In both 
cases we take e = 1.2 x 10 9 m, L = 10 19 m and Te = —s for Ze = 3 x 10 24 m; with these 
values \h\ ~ Jj ~ 10 -15 which is within the expected accuracy of PTA [5]. Similar results are 
obtained for other close values of Te 



or a pulsar nearly aligned (but not quite as otherwise E^fi 1 ^ = 0, although the total timing 
residual is non-vanishing due to the 0(A) time components for TT waves) with the source. 
To understand this enhancement let us analyze the behavior of the integral 

1 = J dx(cos© + sinO), (9) 

with O defined in ([7]) as the pref actors in (|8|) are not relevant for the discussion. The result can 
be expressed as a combination of Fresnel functions, and sines and cosines. In the limit where 
A — > the Fresnel functions go to a constant and the behavior is the usual for trigonometric 
functions. In this respect, the Fresnel functions are responsible for the position and magnitude 
of the enhancement. This is clearly seen when I is plotted^ function of the angular 
separation a between the source and the pulsar. I 2 always shows a maximum, the position 
of which is quite stable under changes of most of the parameters involved. It turns out to 
only depend strongly on the value of A and on the distance to the source. It actually depends 
on the time scales involved rather than on the distance to the source but since the time of 
arrival of the wave to the local system is directly related to the distance, the dependency 
is correlated. This is evidenced in Figure [3] which shows plots of I 2 for different values of 
the frequency, distance to the pulsar, distance to the source and cosmological constant. In 
Figure El a the following reasonable values, Ze = 3 x 10 24 m, w = 10 _8 s _1 , Te = —s and 
L = 10 19 m are used. In b) there is a change in the distance to the pulsar. In c) we change 
the frequency. In d) we keep the distance to the source fixed and use the time at the end of 
an hypothetical 3 year observation. In e) we change the distance to the source one order of 
magnitude (therefore time also changes). Finally in f) the cosmological constant is changed. 
It is clear that the most dramatic changes occur when either the distance to the source or 
the value of the cosmological constant are modified. 

3 We plot I 2 rather than I to deal with a positive quantity. 
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3 Significance of the timing residuals 

Now we would like to make a more detailed study of this possible signal. For that we use the 
ATNF pulsar catalogue [12]. As it is well known pulsars are remarkably stable clocks whose 
periods are known to a very high accuracy, up to 10 -14 s in some cases. However to achieve this 
extreme precision requires some hypothesis that are not appropriate for the physical situation 
we are considering and we will assume the more modest precision of at = 9.6 x 10~ 7 ~ 10 -6 s. 
This value is obtained by averaging the precision achieved for the best measured pulsars 
included in the International Pulsar Timing Array Project, Table 1 in [3]. We are aware that 
only around forty pulsars are monitored with such accuracy for the time being. Any future 
improvements in precision and scope would directly translate into an improvement of the 
results presented in this section. 

For each pulsar we have the galactic latitude ((ft), the galactic longitude (6) and the 
distance (L) to the Earth. We transform these coordinates to (a,/3), where a as already 
explained is the angular separation between the line Earth-GW source and the line Earth- 
pulsar. /3 corresponds to the azimuthal angle of the pulsar referred to the plane perpendicular 
to the line Earth-source. 

The statistical significance of the timing residual will be 



where at is the accuracy with which we are able to measure the pulsar signal period. We 
take at = 10 _6 s as mentioned. The index i running from 1 to N p labels the pulsars included 
in the average. 

In the statistical average we assume an observation time span of approximately 3 years, 
starting at the time the signal is 10 16 s old (time of arrival at our Galaxy). We assume that 
we perform observations every 11 days. That is N t = 101; 10 16 s <T E < 1.00000001 x 10 16 s. 
Since the coalescence times of super massive black hole binaries (SMBHB) can be of the order 
of 10 7 s p] (that is a much shorter time scale than the time of arrival of the perturbation 
to the local system) it is justified to use T E = 2m. . Form Fi gure EH one can also see that 
the position of the enhancement is not significantly altered in the time span of observation. 
While we are aware that 3 years is a short time (most studies consider observational periods 
from 5 to 10 years) we do not intend to present here but a proof of principle and we prefer 
to consider a short period for our numerical analysis. Longer periods of observation will of 
course reinforce the signal. 

We turn to the angular dependence of the significance. In the following a (a) is plotted 
keeping a as a free parameter (note that it is not summed up), that is, using a set of 5 fixed 
pulsars supposed to be exactly at the same angular separation^ from a source the position 
of which we vary < a < tt (the result again is symmetrical for tt < a < 2ir). This could 
be done for any set of five pulsars, since, as shown in the previous section, the position of 
the peak does not depend on the values Lj and However, we used the following set of 
real pulsars which are all close to each other at a distance L ~ 10 20 m. It must be borne in 
mind that although there are over 600 pulsars, making it easy to find clusters with a similar a 
(albeit possibly with very different values of L and j3) , the precision with which they are timed 

4 The chosen pulsars belong to a globular cluster that in principle can not be currently timed with the 
assumed accuracy due to the internal accelerations within the cluster, in this respect this is still a theoretical 
exercise. We thank the referee for pointing this to us. 




(10) 



6 



Pulsars from the ATNF Catalogue 
J0024-7204E 
J0024-7204D 
J0024-7204M 
J0024-7204G 
J0024-7204I 

Table 1: List of pulsars whose Li and j3% we used to calculate cr{a) for an hypothetical source 
at angular separation a. 



can vary widely. The overall magnitude of the significance depends directly on the precision 
of the measured period as well as on the amplitude of the wave. The results presented in 
the following may not be entirely realistic due to the uncertainties in these values but the 
general features of the analysis would remain unchanged if more realistic values (eventually 
available) are used. 



cr (a) 



\ 



1 ^/ g(^,L i ,q,A,3x 10^,10-8,1.2 x lQ9,lQ-35) \ 

ioi 2^Z> I io-6 J • 

i=l j=l \ I 



Length units are given in meters, frequencies in s . We observe a huge peak at a ~ 0.19rad 
(see Figure H|) . If a source is located at such angular separatiorH from the average angular 
position of the five pulsars chosen for observation, the significance could be boosted some 
50 times. Let us compare it to the same calculation taking A = and redshifted frequency 
w e ff = 7i^iy; z ~ 0.008, which is the corresponding redshift for an object 10 24 m away 
calculated using both matter and energy densities. No peak is observed. 

Now we take a list of observed pulsars well distributed in the galaxy. The angles (a,/3) 
are calculated for all of them considering two hypothetical sources of GW. One located at 
galactic coordinates Ogi = 300°, <f>si = —35° and another located at 6s2 = 4°, <p§2 = 10°. 
We order the pulsars from the lowest a to the largest. We group them in sets of five pulsars. 
We consider 27 sets of 5 pulsars; that is a list of 135 pulsars. For each set we calculate the 
significance 



1 Y^j^ HjT^^^a^^Sx 1024,10-8,1.2 x lQ9,lQ-35) \ ^ 



and plot it as a function of the average angle of the set, a& = Y2i=i it with 1 < k < 27. 
Note this is different from pip; here we choose two hypothetical fixed sources and a long list 
of pulsars grouped by their angular separation a to these sources. This could be a realistic 
calculation once real sources are considered. 



5 The angular position of the enhancement corresponds approximately to a stationary point of the phase 
of the wave, i.e. the path along which the phase of the gravitational wave is practically constant, making the 
integrated timing residual maximal. If one considers the wave front of a Minkowskian spherical wave, there 
is no path parameterized as the argument in (4), corresponding to a straight line in space, with a constant 
phase (other that paths pointing to the source, for which the amplitude is zero). Instead, if one considers a 
wave with a frequency and wave vector depending on the space-time coordinates, it is in principle possible 
to find a path where the changes in frequency and wave vector compensate the change in phase, making the 
overall phase constant along that path. It can be seen by deriving the phase of the trigonometric function 
in hij with respect to x and imposing the stationary phase condition that one gets a solution for a/0 that 
depends very weakly on x itself. 
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The results obtained are plotted in Figure [5j In both cases a very noticeable peak is 
observed at the expected angle. 

The reason why the peak for Source 2 is lower than the peak for Source 1 is that Source 
2 is located close but not at the precise angular separation of a real cluster of pulsars. This 
is meant to illustrate that even in that case a significant enhancement of the signal can be 
achieved. 

Finally, the dependency of a on the frequency 

H{T^,L iiai ,^ x 10 24 , w ,1.2 x 1Q9, 10-35) y 

10=6 J > W 

has also been investigated. Some of our preliminary checks indicated that no differences at all 
were observed in the power spectrum when the value of A was changed and that, as expected 
[2J [31 [3], the signal follows a power law a ~ ^. However, let us take a closer look at the 
dependency on the frequency for a short list of pulsars located at the right angular separation 
to observe the peak. We have already seen the significance grows notoriously in this angular 
region. Figure [6] (middle) shows the frequency dependence of the signal for fifteen pulsars 
at the right spot with respect to Source 1. As it can be clearly seen, the signal significance 
grows enormously again for A = 10 -35 s -2 and apparently does not follow a power law. For 
the same short list of pulsars and for A = the signal falls back to smaller values and its 
envelope shapes towards a power law. In Figure (top) we also present the same plot for 
fifteen pulsars located at an angular separation of around a ~ l.lrad, that is away from 
the peak. In this case we see no differences between the different values of the cosmological 
constant as well as a clear power law behavior. The magnitude of the signal is compatible 
with that of the fifteen pulsars at the peak separation when A = 0. 



a\w) 



\ 



iv p ioi 



N„ 



101 

i=l 7=1 



4 Measuring the cosmological constant 

We have seen in the previous that there is an enhancement in the timing residual for a 
particular value of the angle a when GW propagating in de Sitter space-time are measured. 
Among all the dependencies, and when the distance to the source is well known, the most 
relevant appears to be the one related to the value of the cosmological constant A. The 
position of the peak depends strongly on the value of A. It moves towards the central values 
of the angle for larger values of the cosmological constant. 

The values of A as a function of the position at which the peak would be found are plotted 
in Figure [7] (dots) using the positions found in the plots for a (a) (fTTj) for different values of 
the cosmological constant. This calculation was carried out using two independent numerical 
methods in order to make sure that one is free of numerical instabilities (this is a necessary 
precaution as large numbers are involved). 

We argued in Section 2 that the position of the peak is determined by the Fresnel func- 
tions one obtains when calculating the timing residuals. Indeed the integral I in Q, which 
captures the crucial effect, gives a prefactor times a combination of Fresnel functions times 
a combination of ordinary trigonometric functions. The latter are featureless; however the 
prefactor becomes quite large for a specific value of the parameters involved. This particular 
value renders the Fresnel function close to zero and the product is a number close to 2. Away 
from this point the net result is small. 

Using the series expansion of the Fresnel functions at first order we are able to obtain 
an approximate analytical expression for the relation A(a); that is for the value of the cos- 
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mological constant that (all other parameters being fixed) gives a strong enhancement of the 
significance a at a given angle a 



k{a) 



12c 2 sin 4 (f) 



12c 2 sin 4 (f) 



(14) 



({cT E - Ze) cos a + Ze) 



2 — 



which is also shown in Figure [7] (line). We have used the fact that, taking into account the 
duration of a black-hole merger, cTe — Ze- Eq. (|14|) is a clear prediction that could be 
eventually tested. In fact, this effect could also facilitate enormously the detection of GW 
coming from massive binary black holes by carefully selecting and binning groups of pulsars, 
although the possibility of measuring A locally certainly looks to us more exciting. 

Throughout this work we have considered only the effect of A on GW and the way they 
affect pulsar timing residuals and we have neglected the effect of matter or matter density. 
In fact, the main effect of the latter would be through the familiar redshift in the frequency 
of GW. The overall frequency value does not play a crucial role in the previous discussion 
provided that is low enough to be detectable in PTA. It is its dependency on the space-time 
coordinates thats brings about new effects. It is probably useful to remind the reader once 
more that A is assumed to be an intrinsic property of space-time, present to all scales, except 
close to the Galaxy. It would be easy to implement more realistic models in our study, if 
reasonably well-defined ones were available. In fact, these uncertainties constitute strong 
reasons to try to measure A locally. 

5 Summary 

The purpose of this work was to investigate the local effects of the cosmological constant 
for the detection of gravitational waves in PTA. The gravitational wave function is usually 
modeled as a massless, either plane or spherical, wave traveling in flat space-time. The 
expansion of the universe is accounted for by including a redshift in the frequency. Major 
problems are related to modeling the source and assessing the strain of the amplitudes of the 
waves. Here we obviate these by just assuming a spherical wave and focus in the fact that 
the waves propagate in a de Sitter space-time rather than in flat space-time. 

We use a wave solution previously derived in FRW coordinates, which we expect to be 
more realistic than the redshifted usual waves. With this, we calculate the timing residuals 
induced in the signal of known pulsars in our Galaxy, predicting a particular value of the 
angle subtended between the source and the pulsar where an enhanced significance of the 
timing residual is observed. We argue that the position of this peak depends strongly on the 
value of the cosmological constant. This peak is absent when the calculations are carried 
out with usual, Minkowski solution, redshifted waves. We propose two hypothetical sources 
at two distinct positions for which we calculate the timing residuals significance using a real 
set of pulsars. The peak is observed at the predicted angular position. Finally we obtain 
the angular dependency of the value of the cosmological constant using the position of the 
peak for different values of A and analytically from the Fresnel functions involved in the 
calculation. This method could represent an independent way to determine the value of the 
cosmological constant. 

We should end up with a disclaimer. The results presented in this article are by all 
means preliminary. We have addressed the somewhat academic study of an isolated point- 
like source of GW and proceeded to analyze its influence on PTA. In fact, we should expect 
a complete background of sources. Preliminary studies indicate that single source detection 
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may be feasible, but at values of the red shift larger that those considered here [3] . Of course 
these studies do not consider the effects of A discussed in this paper. In addition a full error 
analysis including all unknowns in the appropriate covariance matrix should be performed 
before drawing conclusions on the statistical significance of the effect. Some of these studies 
are only feasible by the PTA collaborations themselves whose interest on the effect presented 
here we hope to have aroused. 
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Figure 3: Integral I 2 plotted for different values of the parameters involved, a) Corresponds 
to the reasonable values w = 10~ 8 s -1 , L = 10 19 m, Ze = 3 x 10 24 m, A = 10 _35 s -2 and 
Te = 10 16 s. b) Change in pulsar distance to L = 10 21 m. c) Change in frequency to 
w = 10 s . d) Change in time to Te = (10 16 + 10 8 )s. e) Change in time and distance to 
the source to Ze = 3 x 10 23 m and Te = 10 15 s. f) Change in the cosmological constant to 
A = 10" 36 S - 2 . 
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Figure 5: Plot of a k (a k ), k = 1,27. A = 1(T 35 
squares to Source 2. Full range is showed on top, zoom on the lower values for A = 10 -35 
and comparison to A = show on middle and bottom respectively. 
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Figure 6: a{w) for 15 pulsars away from the peak angular region for Source 1 (top), the solid 
line corresponds to A = 10 _35 s -2 and dots correspond to A = 0. a{w) for 15 pulsars at the 
peak angular region for the same source (middle). Solid line corresponds to A = 10 _35 s -2 
and the data close to the horizontal axis corresponds to A = 0. Zoom on the A = case 
(bottom). 
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Figure 7: A(a) obtained numerically from the positions of the peaks in the a{a) plots for 
different values of the cosmological constant (dots) and obtained analytically from an ap- 
proximation of the Fresnel functions involved in the timing residual (line). 
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